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Theoretical calculations of one- and two-hole spectral functions for the 16 nu- 
cleus are still failing to describe some of the important features observed exper- 
imentally. Of critical importance for the solution of these issues is to obtain an 
■ appropriate description of the interplay between hole-hole and particle-hole ex- 

' citations. A formalism is reviewed here that allows the consistent treatment of 

| such phonons at a Random Phase approximation level. Although the application 

of this formalism to 16 is still in the implementation stage, some preliminary 
results are discussed here. 
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INTRODUCTION 



As a result of correlations between nucleons inside the nuclear medium, the occupation probabil- 
ity of single-particle shells are lowered with respect to the mean-field prediction. Experimentally, 
this can be observed as a quenching of the absolute spectroscopic factors for the knockout of a 
CN ■ nucleon from a given shell. Studies of (e, e'p) reactions have allowed the determination of spectro- 

scopic factors in many closed-shell nuclei Jl|, ||, [| demonstrating that the removal probability for 
nucleons from these systems is reduced by about 35% with respect to the shell-model predictions. 

For the 16 nucleus, the experimental spectroscopic strength Q for the knockout of a proton 
from both the pjy 2 and p 3 / 2 shells corresponds to about 60%. The inclusion of relativistic effects 
in the analysis of data is espected to change these results by not more than a few % Theo- 
retical calculations for this nucleus seem still far from being able to reproduce quantitatively the 
experiments and give numbers that depend on the approximation scheme employed. Variational 
calculations, focussed on short-range correlations, yield about a 10% of strength removal [0, 0, || 
before the effects of the center of mass motion are included. Center of mass corrections are known 
to raise the spectroscopic factor by about 7% |j, thus worsening the agreement with data. An 
improvement of these results seems to require a treatment of low-energy correlations with the inclu- 
sion two-hole-one-particle (2hlp) states [1C, 11]. This would be in line with the Green's function 
calculations of Ref. | |12[ |, in which 2hlp states were taken into account and a reduction of about 25% 
was obtained, still neglecting center of mass corrections. A comparison with the experimental result 
(and their uncertainties) suggests that about 20% of the observed single-particle strength removal 
remains unexplained. Nevertheless, it is important to note that the results of Ref. [12] indicate 
the low-energy correlations as an essential ingredient, needed for a complete understating of this 
puzzle. 

The merit of the calculation of Ref. [12] lies in the fact that summing 2hlp states and their 
interactions to all orders, one can achieve a simultaneous description of the effects of both particle- 
hole (ph) and hole- hole (hh) collective excitations, including the interplay between them. These 
effects, though, were accounted for only at a Tamm-Dancoff approximation (TDA) level. In order to 
account for the coupling to collective excitations that are actually observed in 16 O it is necessary to 
at least consider a Random-Phase approximation (RPA) description of the isoscalar negative parity 
states [O]. To account for the low-lying isoscalar positive parity states an even more complicated 
treatment will be required. Sizable collective effects are also present in the particle-particle (pp) 
and hh excitations involving tensor correlations for isoscalar and pair correlations for isovector 
states. 

Taking a look at the experimental spectral function p], one can also notice that the total strength 
for proton emission of a p%/2 proton is experimentally aistributed over a main peak and two small 
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fragments at slightly higher missing energy. This splitting is not described by the calculations of 
Ref. [12] where only one pole is generated that contains most of the single-particle strength. 

The results of Ref. (T^| have been used as a starting point to study the 16 0(e,e'pp) reaction [14, 
(E5[ |. In these works, the two-hole spectral function was obtained by coupling two dressed single- 
particle (sp) propagators and then allowing them to interact to all orders through a ladder type 
equation. The distortaion of the two-hole overlap function caused by short-range effects was also 
included by adding defect functions associated with a G-matrix interaction [16| which was computed 
for the finite 16 nucleus. These calculations led to a proper description of the experimental cross 
section for two proton emission [ ]l7| , |i~8| ] and yield amplitudes for the transition to the ground state 
of 14 C and to some of its lowest excited states. 

Nevertheless a completely satisfactory description is missing here also. Recent developments of 
NN knockout reactions allow to distinguish several states of the final nucleus and in particular it 
is now possible to disentangle the two 2 + excited states of 14 C at 7.0 and 8.32 MeV pcfl . In 
spite of this, the theoretical calculations of Ref. [14] can generate only one 2 + excited state. It is 
interesting to note that the transition to both of the 2 + states can be interpreted as the knockout 
of two protons from a p±/ 2 and a p 3 / 2 state. Although this has not been directly investigated yet, 

it comes natural to suppose that the fragmentation of the p 3 / 2 peak seen in the 16 0(e,e'p) reaction 

could play some role in the generation of the 2 + doublet in 14 C. Moreover, it must be noted that the 
strength of the pi/ 2 and p%/ 2 peaks, generated in the calculations of the one-hole spectral function, 
enters directly in the calculation of the two-hole amplitude. 

It is therefore evident that one- and two-hole spectral functions are closely related, the correct 
description of the first being relevant for a successful description of the second and vice versa. As 
a consequence, a formalism in which both quantities are generated together in a consistent way 
is highly desirable. This can be achieved in the the framework of self-consistent Green's function 
theory (SCGF) in which the equations of motion are expressed in terms of fragmented propagators 
and are iterated until self-consistency is reached. 

At the same time, it is also desirable to extend the calculations of Ref. [12] to include the 
RPA description of pp (hh) and ph phonons. This is not completely trivial. For example, the 
naive summation of diagrams containing both pp and ph phonons leads to serious inconsistencies. 
This approximation is depicted in Fig. [j]. The last of the three diagrams on the right hand side 
is already contained in each of the other two and must therefore be subtracted to avoid double 
counting. This subtraction introduces spurious poles in the Lehmann representation of the self- 
energy and generates meaningless (not even normalizable) solutions of the Dyson equation. In 
addition, each of the first two terms in Fig. [I] ignores the Pauli correlations between the freely 
propagating line and the quasi-particles forming the phonons, as noted in |21[. 

To achieve the inclusion of RPA effects, one needs to pursue a formalism in which the contri- 
bution of the pp and ph phonons to the self-energy are summed to all orders, like the example of 
Fig. [|. This, woul d av oid the subtraction of the second-order diagram. Such a formalism has been 
proposed recently [22| with the intent of attacking the above issues. In this approximation scheme 
the effects of both pp (hh) and ph phonons can be evaluated at an RPA level and then summed 



up to all orders by means of the Faddeev technique |23, 24]. The treatment of Pauli correlations 
is also improved over methods that include ph RPA phonons directly in the self-energy, since all 
exchange terms at the 2plh level are consistently included. 




FIG. 1: Example of an approximation for the self-energy. Although this approximation contains both ph 
and pp correlations it would generate incorrect results, due to the need for subtracting the 2nd order term 
to avoid double counting. 
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An issue that one needs to be aware of in the application of the formalism of Faddeev equations is 
the apparence of spurious solutions. The Faddeev approach consists in substituting the Schodinger 
equations with a set of three coupled differential equations. As a result, more solutions are generated 
and about 2/3 of them are spurious. In the general three-body problem, though, the properties of 
such spurious states are well known and they can be eliminated efficiently 25]. The situation is a 
little more complex if the Faddeev formalism is applied in a many-body context, with the intent of 
studying the motion of quasi-particle and quasi-hole excitations. Here, in particular, the fulfillment 
of closure relations for pp and ph amplitudes is related to the behavior of the spurious Faddeev 
eigenstates p^] . This makes the choice of the interaction boxes to be used a critical problem, since 
without a proper treatment of this relation the spurious solutions will irremediably mix with the 
physically meaningful ones. 




FIG. 2: The application of the Faddeev 
formalism allows to sum the effects of pp 
(hh) and ph RPA phonons to all orders, 
generating diagrams like this one. As a re- 
sult of the all-orders summation, no incon- 
sistent poles will appear in the self-energy. 



The Faddeev approach represents a step further with respect to the calculations of Refs. [12|, [Tij 
and it is our aim to apply it to the study of one- and two- hole spectral functions of 16 O. In this 
contribution, the basic formalism is briefly reviewed. Some details and some technical issues (as 
the elimination of spurious solutions) are described more completely in Ref. 22] and will only be 
mentioned. The possible implementation of the SCGF scheme is also a key feature of the present 
formalism and it will also be reviewed in some detail at the end of Sec. |2 C| 

The actual application to the nucleus of 16 O is currently under way but has not been completed 
and will be presented in a future paper. Nevertheless, some preliminary results have been generated 
for a simplified model and they will be described in the last part, as an example. 



2. FADDEEV FORMALISM FOR 2H1P AND 2P1H MOTION 



The relevant quantities for the study of one- and two-nucleons knockout reactions are the one- 
and two-hole spectral functions 

s*h = EK**" 1 l^l*o>l a % + ^*" 1 -^). W 

k 



Sapi") = J2\(*i~ 2 \c a c^)\ 2 5(co + Et 2 -E£). (2) 

k 
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These can be obtained from the diagonal elements of the one- and two- body Green's function 
respectively, which in their Lehmann representation take the form [26, 27] 



V w - eZ + zry Y w ~ £ k ~ l V 

and 

n+ u - e n X + til k- u - £ k- ~ l V 

In Eq. (|), X£ = (#£ +1 [4 ) (y% = (^ _1 |c a |*^)) are the spectroscopic amplitudes for the 
excited states of a system with A + 1 (A — 1) particles and the poles e+ = E^ +l — Eq {e7 = Eq — 
E^ -1 ) correspond to the excitation energies with respect to the ^4-body ground state. Completely 
analogous definitions are employed in Eq. (|j) for the transition amplitudes and excitation energies 
related to the states of a system with A + 2 {A — 2) particles. 

In the nuclear case, a strong coupling exists between the sp degrees of freedom and both collective 
low-lying states as well as high-lying states. The latter coupling is related to the strong short-range 
repulsion in the nuclear force. The resulting fragmentation of the sp strength (as observed in 
experimental data) suggests that this feature must already be included in the description of these 
couplings. This should be done by expressing the equations used to evaluate the propagators g a p{oj) 
and gjwapiw) in terms of the solution g a p{oj) itself. This self-consistency feature also emerges in 
an exact formulation, involving the coupling to two-, three- and A-body propagators, which can be 
derived using the equation of motion method |2^] . 

Keeping in mind this approach, we compute g a p{^) as a solution of the Dyson equation 

M") = 9 ( $M + , (5) 

where 2*^(0;) is the irreducible self-energy. Here and in the following, sums over repeated indices 
are implied. 

By considering the equation of motion for g a p(u)), one obtains that S*^(cj) can be written as 
the sum of two terms 

= S a/f + jVuX^v RpvX^Seiu) VyS,j3e , (6) 

where represents the Hartree-Fock part of the self-energy, which can be computed straightfor- 
wardly from the solution g a p(uj) itself. The V a \^ u represent the antisymmetrized matrix elements 
of the residual interaction. 

The non trivial element to be inserted in the calculation is the propagator R^vX^Sei 1 ^) , appearing 
in the last term of Eq. (JBJ), which describe the motion of excitations consisting of at least two- 
particle-one-hole (2plh) or 2hlp states. It contains the sum of all six points diagrams that are 
one-particle irreducible, i.e. that cannot be separated by cutting a single line. 

In general, R^uX^Sei^) is the solution of a Bethe-Salpeter type equation involving propagators 
depending on four times (i.e. three energies) as well as four- and six- point kernels p9p . This is 
appears too complex to be solved numerically and one needs to find an appropriate approximation 



for the propagator R^uX^Sei^)- The formalism proposed in Ref. 2^1 aims to compute this object 
keeping the effects of RPA phonons but at the same time restricting oneself to equations involving 
only two-times propagators, thus producing a set of equations that one can deal with. This is 
achieved by applying the Faddeev equations technique. 

Also, the resulting formalism splits up in two separate expansions for the 2plh and the 2hlp 
components. Although the hole spectral function is of primary interest for comparison with exper- 
imental data, it must be stressed that both 2plh and 2hlp components are needed to generate the 
self-consistent solution for the sp propagator. Since the formalism involved is the same for both 
components, we will describe only the forward-going (2plh) expansion in the rest of this section. 
The equations for the 2hlp case are completely analogous. 
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A Faddeev equations 

As a convention, we assume that the first two lines propagating through the expansion of 
R^uX^Sei 10 ) represent particle excitations while the last represents a hole, as in Fig. H Following 



standard notation in the literature [3C], the i-th Faddeev component R^l^ aB^y wm re P resen t the sum 



^ii/X,af3 r y 

of all diagrams ending with an interaction between legs j and k, with (i,j,k) cyclic permutations 
of (1,2,3). Obviously, the component R^ 3 ' refer to the interactions in the pp channel, while the 

two components i£W and Ry) will refer to the interaction between a particle and a hole. The last 
two are trivially related to each other by the exchange of the first two indices. 

Also, for reasons to be mentioned in the next subsection, the By) need to be redefined in such 
a way that their matrix elements also depend on the indices (n, n', k), which label the fragments 
of the external propagator lines (compare with the definitions below Eq. (|3|) ). This implies that 
the eigenvalue equations will involve summations on both the sp indices (a, (3, 7) and the ones 
corresponding to the fragmentation, (n a , n^, /c 7 ). The 2plh propagator and its Faddeev components 
are recovered only at the end of calculations by summing the solutions over all values of (n a , rip, 
fc 7 ) and (n M , n u , k\). 

Putting together all the above ingredients, the resulting approximation to the Faddeev equations 
can be rewritten in a way where all the propagators involved depend only on one energy variable 
(or two time variables). The 2plh part of this expansion can be written as follows 

o(«) ( \ 

fi,n^vn I/ Xk\,ctn a l3np'Yk~ l \ I 



fj,n ll i/n u Xkx,an a Pni3'Yk 1 \^) un v ^in li Xkx,an a l3np'yk~ t y^) 



+ C°i/n,^n M Afc A ,/i'n^'<A'fc^ ( w ) n' vi i' n'^ X 1 k' x ^" n^v" n'lX" k'^ ( W ) 

X (^"ri^''<A''fc^an Q /3n /3 7fc 7 ( a; ) + R f Jn'/ t v''n'iX''k'l,an a Pn ^k^ U ') > «' = 1> 2, 3 . (7) 

In Eq. (0), G° is the forward-going part of the 2plh propagator for three dressed but noninteracting 
lines. Using the notations introduced after Eq. (||) we have 

G im li vn v \kx,a.n a PnB'yk.y( u) ) = <W,™/3 $k\fa — , + + n . • ( 8 ) 

<JJ — [tn a T £n — £ fc7 J -r IT/ 

The Faddeev vertices T„)„, „,„, ,,,, „//_//, ,„„,, \nurr(w) contain the collective excitations in the pp 

V Tl u \l Tl^A K^ifl TV^y 'V A 

and ph channels and can be evaluated in RPA, as explained in the next subsection. 

Once the Faddeev equations are solved and the usual representation of the Faddeev components 

R^lx a/3j * s recovered by summing over the indices (n a , np, fcy) and (n^, n v , k\), the complete 2hlp 
propagator can be obtained as the sum of all the components: 

Rnv\,a/3j(v) = R ^uX,aP^ u ) ~ (^AW'A.a/^yM ~ G^a^M) , ( 9 ) 

i=l,2,3 

The free forward-going propag ator G° was already included in the Eqs. (|7|) and therefore is to be 
subtracted in Eq. (g) to avoid double countings. It must be stressed though that this subtraction 
does not generate any inconsistencies in the self-energy, as it was the case for the approximation of 
Fig. |. Rather, this subtraction is required to exactly cancel similar, troublesome, terms that are 
included in the single Faddeev components. 

B Faddeev vertices 

The interaction matrices ^j^/ a p(ui) for the pp channel and n^f^(a;) for the ph channel, are 
obtained by solving the ladder equation in the Dressed RPA (DRPA) scheme, which employs dressed 
propagators as they are obtained by solving the Dyson equation (ph. 



© submitted to 5 th Workshop on "c-m induced Two-Hadron Emission", Lund 2001. February 8, 2008 - 6 



For the pp case, the DRPA equation is graphically depicted in Fig. 0. The gamma matrix 

( \ 

r^f^(oj) obtained in this way is the generalization of the two-body T-matrix for two nucleons and 

represents our approximation for the two-particle propagator of Eq. (Q) , when the external legs are 

amputated. Completely analogous considerations apply to the ph interaction matrix n^L^a;). 

The RPA phonons generated in this way depend on dressed propagators and can be included in 
the Faddeev formalism outlined above. It is worth noting that discarding the last diagram on the 
right hand side of Fig. B and in the corresponding ph counterpart, the whole formalism reduces to 
the TDA expansion of Ref. [12]. 



p(pp) 



■• + 



p(pp) 



p(pp) 




FIG. 3: 

tion for 
trix. 



DRPA equa- 
the Y^v) 



ma- 



The vertices to be employed in the Faddeev equations (0) will consists of the interaction boxes 
just described to which a freely propagating dressed line is added, as displayed in Fig. ||. The exact 
formula for the pp channel is as follows: 



P 



(3) 

fin ll i/n^Xkx,an a f3nfj'yk-y 



To 



v, 



5 W 



fJ,U 



A 



a/3 



+ E 

k- 



£ k X )+ lT l 



e nij + e k x + £ k- 



(4: 



et)(4Z 



(10) 



were the A n+ (A fc ) amplitudes are computed through the DRPA equation of Fig. |3| and are related 



to the amplitudes of Eq. @ by A™+ = V ma p (A*~ 

The complication of defining the interaction vertex not only in terms of the sp indices but 
also in terms of the respective fragmentations (n M , n u , k\) and (n a , np, k 7 ) is a consequence of the 

two main requirements imposed: that is, has to depend on only two times and the free line 
added in Fig. |4]must be a dressed propagator. As it can be seen from Eq. (10), the energy poles , 
relative to the fragments of this line, appear in the equation mixed with the corresponding poles of 

i(pp) 



the r^a j a(w)'s phonons. This situation couldn't be handled without the above prescription |22j, |35|j. 
Instead, with this procedure it is still possible to write an expression for the diagrams of Fig. || in 
the form G 0> (uj) r^ 3 ^(a;) G 0> (to) or, equivalently, to write the Faddeev equations (f7|) . 

As a last remark, we note that the last term on the right hand side of Eq. (JlQ) comes from 
the last diagram of Fig. H. This term contains large energy denominators and it is expected to 
give relatively small contributions. Nevertheless, it can be proved that it is essential to include it 
in order to impose the appropriate behavior of the spurious solutions generated by the Faddeev 
equations (and eventually to discard them) |22|. Thus also this diagram needs to be kept. 



C Self-consistent Scheme 



The set of Faddeev equations (0), together with the expressions (|5|) and (q), furnish a scheme 
in which one- and two-body spectral functions can be evaluated in terms of an already dressed 
single-particle Green's function. This allows the possibility of including, already from the start, 
the features of the single-particle motion that go beyond a mean-field description. Once such 
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FIG. 4: Here Ar> and AT < are the forward- and backward-going part of the energy dependent contribution 
to the pp DRPA vertex (0). The contribution of these three diagrams can be factorized in an expression 

of the form G 0> G 0> only after having redefined the propagators G 0> and I^ 3 -* to depend also on 
the particle and hole fragmentation indices (n,n' ,k). The last diagram has a smaller effect on the physical 
solutions of the problem, although it is essential for the elimination of spurious solutions. 



features are included in the input propagator, their effects on the one- and two-body motion will 
be automatically included in the calculations. It is important to observe that the expansion of the 
irreducible propagator R^x^Ssiuj) in Eq. @, derived employing the equation of motion, is given in 
terms of the real one-body propagator. Thus, in principle, it is the exact solution of the many-body 
problem that is supposed to be employed in the Faddeev calculations described above. 

The approach of self-consistent Green's function theory consists in starting with an approxi- 
mation for the input Green's function (usually an Hartee-Fock propagator or other independent 
particle model). From this, the pp (hh), the ph and the Faddeev propagators can be evaluated with 
the formalism described. Then, the solution of the Dyson Eq. (||) will give a better approximation 
of the s.p. propagator, that can be employed in the second calculation. As Shown in Fig. S the 
whole procedure is iterated many times until consistency is found between two successive solutions 
(i.e. between the input propagator and the solution that it generates). 

The attractive feature of this procedure is not only restricted to the fact that the effects of 
fragmentation are included in the calculations. But also, expressions for the one- and two-hole 
motion and the ph phonons (which give information on the excitation spectrum of the system) 
are generated all at the same time, the relative effects of each one on the others being taken into 
account. 

Obviously, the exact form of the single-particle propagator is very complex, containing many 




FIG. 5: Iteration scheme used to reach self-consistency between the input propagator and the result of the 
Dyson equation. At each iteration new approximations for the one- and two- body propagators (|3|) and (jj) 
are obtained. 
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TABLE I: Single particle energies e° a and interaction strengths used in Eq. ([Tl|). The sp energies were 
computeded from Ref. fll2| , while g Q and g v were fitted to roughly reproduce the excitation energy of the 
3 _ ph phonon. 



a ls 1/2 lpx/2 lp 3 / 2 14/2 2s 1/2 


14/2 1/7/2 


2p 3 /2 


I/5/2 


2^1/2 


159/2 


e° a (MeV) -35.0 -18.2 -9.3 -1.8 -0.1 


4.4 17.4 


16.0 


23.5 


17.7 


30.0 


g Q = -l.UMeV/fm 4 




9 V = + 


15.0 MeV 







poles and a continuum. Such richness of details cannot be easily included in the above calculation. 
In our case, though, we are interested in the low energy behavior of the single-particle motion and 
the few poles that contain the main sp strength are expected to be much more relevant than the 
others. One possible prescription to limit the number of poles in the input propagator at each 
iteration consists in keeping only the most significant peaks and collapsing all the others in an 
effective pole. This procedure was already employed in applying the SCGF method to the study 
of pairing in superfluid nuclei |32j . 

The effects coming from short-range correlations, can be taken into account by substituting the 
bare NN potential V a fi^b with a G-matrix interaction. In the same way, the effects of short-range 
interaction can be included in the two-hole spectral function by adding the relative defect functions 
as done in Ref. [0, HI. 

3. QUADRUPOLE-PLUS-PAIRING MODEL 

The implementation of the Faddeev RPA formalism available at the moment allows us to make a 
test with a simplified model for the 16 nucleus. We considered a model space consisting of harmonic 
oscillator wave functions up to four major closed shells (2p and If states) plus a g 9 / 2 single particle 

state. Within this space, we choose a so called "pairing-plus-quadrupole" Hamiltonian |^3[, in 
which one a body potential U is employed to generate the correct single-particle energies and the 
residual interaction in composed of a quadrupole-quadrupole Vq and a pairing Vp term: 

H = U + Vq + V v 

= E £ ° 4^ + g Q ■ Q ] Q ■ + g P P ] P , (11) 

a 

where P = J2 a c a c a and Q = J2 a ,/3 ( a l r ' 2 ^2,m(^r)|/3)c^c i (3 are the pairing and the quadrupole 
operators, the symbol ":" represents the normal ordering and a is the time reversal of state a. An 
oscillator parameter of b = 1.76 fm was used in computing the quadrupole matrix elements. Both 

Vp and Vq are independent of the spin and isospin degrees of freedom. 

The values of e„ for all the sp states were computed in Ref. [12] in a Brueckner-Hartree-Fock 





a 

s close to 



The 4 



calculation, based on a G-matrix interaction |l(| derived from the Bonn-C potential p4 
used here are are taken from that reference but with some slight modification of the leve 
the Fermi energy, in order to roughly reproduce the experimental missing energies in the Faddeev 
RPA calculation. The exact values employed here are given in table |. 

The remaining two parameters, g Q and g P , were fitted to other experimental energies. In 

particular, g Q was chosen to roughly reproduce the 3~ phonon of 16 0, while g P was fitted aiming 
to reproduce the experimental missing energy of the 14 C ground state. The latter was found to be 
not very sensitive and was then fixed to the value given in table |. 



Using the Hamiltonian (11), the self-energy was computed employing the Faddeev formalism 
with pp (hh) and ph phonons evaluated in both TDA and RPA schemes. As already noticed, the 
TDA approximation was already employed in the calculation of Ref. O] . 
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TABLE II: Results for hole spectroscopic factors of the main p\/2 and P3/2 peaks, obtained with Hamilto- 



nian (11) 





TDA 


RPA 


Pl/2 


.632 


.588 


P3/2 


.668 


.647 



As seen from table [□], the results obtained for the p\i% and p 3 / 2 spectroscopic factors are 4-5 % 
smaller in the RPA case and seem to go in the expected way. 

It's worth to recall that the results from this model have been obtained within a single iteration, 
without the effects of short-range correlations and with a oversimplified interaction. Thus the 
values obtained for the absolute spectroscopic factors are not really meaningful and are expected 
to change significantly in a realistic calculation. Self-consistency can play a relevant role too. 

Nevertheless, a comparison between the results obtained in TDA and RPA calculation seems 
reasonable and can give some idea of the relative importance of RPA effects. 



4. SUMMARY 



The present theoretical description of the distribution of spectroscopic strength at low energies 
lacks important ingredients for a successful comparison with experimental data. One of these 
ingredients is a proper description of the coupling of sp motion to low-lying collective modes that 
are present in the system. Recent calculations for 16 [12], for example, only include a TDA 
description of the se c ollective modes. This contribution outlines a new method that has been 
proposed recently p2j that aims to study the influence of pp and ph RPA correlations on the sp 
propagator for a sysiem with a finite number of fermions. This method is formulated in the context 
of SCGF theory by evaluating the nucleon self-energy in terms of the 2plh and 2hlp propagators. 
The description of the 2plh (or 2hlp) excitations has been studied by using the Faddeev formalism, 
which is usually applied to solve the three-body problem. As mentioned in the text, the detailed 
choice of the form of the interaction matrices as well as the formulation of the problem are imposed 
by practical issues, like the existence of spurious solutions. 

It should be noted that the Faddeev approach outlined here, has also the advantage that it can 
be naturally extended to the inclusion of more complicated excitations like the extended DRPA |35|] . 

The application of this formalism is underway for the nucleus of 16 O. Although no physically 
meaningful results have been obtained yet, the Faddeev RPA equations have been applied to a 
simple model, demonstrating the feasibility of these calculations. The results obtained with this 
model yield somewhat stronger quenching of the py 2 and P3/2 spectroscopic factors when RPA 
correlations are included, thus going in the expected direction. 

A realistic self-consistent calculation in which a G-matrix derived from realistic interaction is 
employed is underway and the results will be reported elsewhere. 
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